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We discuss the extraction of Os using isovector hadronic r decay data and sum rules constructed specifically to 
suppress contributions associated with poorly known higher dimension condensates. We show, first, that problems 
with the treatment of such contributions affect earlier related analyses and, second, that these problems can be 
brought under good theoretical control through the use of an alternate analysis strategy. Our results, run up to 
the Uf — 5 regime, correspond to Qs(M|) = 0.1187 ± 0.0016, in excellent agreement with the recently updated 
global fit to electroweak data at the Z scale and other high-scale direct determinations. 



1. Introduction and Background 

The strong coupling, as, at some convention- 
ally chosen reference scale, is one of the funda- 
mental parameters of the Standard Model (SM). 
Its value, in the Uf = 3 regime, can be extracted 
using hadronic r decay data as a consequence of 
the finite energy sum rule (FESR) relation 

/ w{s)p{s)ds — — (f w{s)II{s) ds {1) 

Jo 27ri Jis\=so 

which is valid for any analytic weight, w{s), and 
any correlator n(s) without kinematic singulari- 
ties. In Eq. IH), p{s) is the spectral function of 
n(s). The basic idea is to use experimental spec- 
tral data on the LHS and, for sufficiently large 
So, the OPE representation of 11 (which involves 
as) on the RHS. The region of applicability of the 
OPE is extended to lower sq when w{s) satisfies 
the condition w{s — sq) — 0, which suppresses 
contributions on the RHS from the region of the 
contour near the timelike real axis 1 2 . 

Experimental input for the LHS of Eq. ^ is 
available because, in the SM, the kinematics of 
T decay allows the inclusive rate for decays me- 
diated by the flavor ij — ud, us, vector (V) or 
axial vector (A) hadronic currents to be written 
as kinematically weighted integrals over the spec- 
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tral functions Pyjj^.j^jis) of the spin J = 0, 1 com- 
ponents of the relevant current-current two-point 
functions [3]. Explicitly, with y-r = s/m^ and 
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with Vij the flavor ij CKM matrix element. 
Sew a short-distance electroweak (EW) correc- 
tion, and p^y'l^ljis) = + Pf/Anji'')- 

For ij — ud, apart from the tt contribution 
to PaL' PyL(s) and p^Aludi^) are numerically 
negligible, being proportional to 0{[md T 
The sum of flavor ud V and A spectral functions, 

PvX^A-udi^)^ *^an thus be extracted from the dif- 
ferential decay distribution dRy+A-ud/ds, for all 
s < ~ 3.16 GeV'^. Further separation into V 
and A components is unambiguous for mr states, 
but requires additional input for KKwr in > 0) 
states, making errors on the experimental distri- 
bution smallest for the V-l-A sum. 

Given the spectral functions Py^^ v+A-udi 
FESRs for the related correlators, H^.^^-*, with 
T = V,A,V + A, are straightforwardly con- 
structable. For the scales Sq > 2 GeV^ considered 
here, the OPE representations of these correlators 
are strongly dominated by the dimension D = 
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contribution, which is entirely determined by as, 
converges well, and is known to 0{a'^) FJl. The 
resulting FESRs are thus well adapted to the de- 
termination of as- To optimize the precision of 
this determination, however, care must be taken 
in evaluating the small, residual higher D con- 
tributions, a ~ 1% determination of as{M'^), for 
example, requiring control of higher D contribu- 
tions at the level of ^ 0.5% of the D — term [5]. 

For ij = ud and sq > 2 GeV^, D — 2 con- 
tributions arc numerically negligible, being ei- 
ther 0(TO?,.d) or 0{a1m1) [S]- £> = 4 contribu- 
tions are, to very good accuracy, determined by 
the RG invariant condensates {mi££), {rrisss) and 
{aG'^), for which phenomenological input exists 
(see Ref. [7] for the explicit forms of these con- 
tributions, and Ref. [51 for details of the conden- 
sate values employed). D > 6 contributions are 
more problematic since the relevant condensates 
are poorly known or phcnomenologically undeter- 
mined. We deal with these contributions by defin- 
ing effective condensate combinations, Cg, Cg, • • •, 

such that [n(g2)]°^f = Ed=6,8,.. ^^/Q"" (up 

to logarithmic corrections) and fitting these quan- 
tities to data. This process is greatly facili- 
tated by working with polynomial weights wl-s) = 
J2m=o bmV"^ defined in terms of the dimensionless 
variable y = s/sq. For such weights, the inte- 
grated D > 6 OPE contributions have the form 
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allowing contributions of different D to be distin- 
guished by their differing so dependences. 

2. Problems With Existing Analyses 

Existing analyses are based on the approach 
pioneered by ALEPH and OPAL [8l9]. In this 
approach, OPE contributions with D > % are 
assumed safely negligible for all weights em- 
ployed, and the quantities as{m'^), (aG'^), Cg 
and Cg are fitted using the sq — values 
of the spectral integrals corresponding to the 
{km) = (00), (10), (11), (12) and (13) "(fcm) spec- 
tral weights", u;('='")(y) = (1 - y)''y"^w^°°\y), 



where w^'^^^y) = (1 — y)^(l -|- 2y) is the kine- 
matic weight occuring on the RHS of Eq. 
ALEPH |8ll0lllj performed this fit indepen- 
dently for each of the V, A and V-f A channels, 
while OPAL [9] performed independent fits for 
the V+A and combined V,A channels. A poten- 
tial problem with the assumption that all I? > 8 
contributions can be safely neglected is the fact 
that w^'''™) has degree 3 + k + m which, from 
Eq. ([3]), implies that contributions with D up 
to 16 are, in principle, present in at least one of 
the FESRs considered in the ALEPH and OPAL 
analyses. Since only the single sq value sq = 
is employed, there is no way to prevent the fit 
from adjusting to the presence of any neglected, 
but non-negligible, D > 8 contributions by shift- 
ing the lower D parameters determined in the fit 
in such a way as to compensate, as best as possi- 
ble, for the missing terms. Such a problem with 
the fit can only be exposed by studying the same, 
or related, FESRs over a range of sq, where the 
different scaling with sq of terms of different D 
will become operative. 

A simple way to test for the presence (or ab- 
sence) of such problems in the ALEPH and OPAL 
fits is to consider the SQ-dependent fit-qualities, 
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where T = V,AotV + A, I^pe{so) and I^p^M 
are the OPE and spectral integrals appearing, 
respectively, on the RHS and LHS of the corre- 
sponding w(s)-weighted FESR, and SI^p^^{so) is 
the error on /]^ec('5o), determined using the ex- 
perimental covariance matrix for dRx-ud/ds. Be- 
cause of strong correlations between values corre- 
sponding to the same w{s) but different sq, a fit- 
ted version of the OPE representation should be 
considered reliable only if |i^^(so)| remains < 1 
for a range of sq < m^. 

In Ref. [S] this condition was shown to be far 
from satisfied for any of the spectral weights em- 
ployed in the ALEPH and OPAL analyses. An 
illustration of the problem is provided in Fig- 
ure [U which shows the Fy (so) corresponding to 
the 2005 ALEPH final data and OPE fit for four 
weights, W2(y) = (1 -y)^, w^iy) = 

1 — |j/ + and w{y) = y{l — y)^, all having 
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degree < 3 (and hence OPE contributions only 
up to D = 8). If the fitted values for the I? < 8 
OPE parameters obtained by ALEPH are reli- 
able, one should find an sq window below hav- 
ing \Fy{so)\ < 1 for all four weights. The results, 
given by the light lines in the figure, show that no 
such window exists. In fact, for the weights W2, 
W3 and j/(l — y)^ not employed in the original fit, 
the fit quality is poor even at sq — rn^. 

The problem seen in the Figure could be due 
either to contamination of the Z? < 8 OPE fit pa- 
rameters by neglected, but non-negligible, Z? > 8 
contributions, or to OPE breakdown. One may 
test the latter possibility by performing alternate 
fits in which potential D > 8 contributions, where 
present, are explicitly taken into account. The 
results of such fits, discussed in the next section, 
yield alternate OPE representations in excellent 
agreement with the corresponding spectral inte- 
grals over a range of sq , both for the weights em- 
ployed in the fits and for related weights with 
OPE representations determined by the same set 
of OPE parameters. The resulting alternate fit 
qualities, Fy (sq), for the four degree < 3 weights 
already discussed above, are shown by the dark 
lines in Figure [1] The results clearly show no ev- 
idence for OPE breakdown. 

3. Alternate FESR Analyses 

In what follows, we employ the updated char- 
monium sum rule determination of (aG^) |12]. 
The gluon condensate term dominates the D = 
4 OPE contribution. Details on the input for 
the small corrections proportional to {meii) and 
(tosSs) may be found in Ref. [5]. 

The D = contributions are evaluated using 
the expression for the D — Adler function series 
from Ref. [4 . An 0{a^) contribution, employing 
the estimated value for the corresponding coeffi- 
cient from Ref. j4j, is included for our central fit. 
We consider both the contour improved (CIPT) 
and fixed order (FOPT) determinations of the in- 
tegrated D = sum. The reference scale Uf = 3 
coupling needed in the evaluation of the CIPT 
and FOPT sums (taken to be as{rn'^)) is a pa- 
rameter to be determined in the fit. 

Since, for V and A correlators, OPE breakdown 



is expected (and observed) to set in for sq below 
2 GeV^ |13|14j . a limited window of sq values 
is available for use in fitting and the unknown 
Cd>4- It is thus convenient to work with FESRs 
based on the weights 

-.(y) = l-^. + ^.^ (5) 

which, like w^^°\y), have a double zero at s = sq 
(y = 1). From Eq. ^ we see that the wn 
FESR involves only a single integrated unknown 

D > 4 OPE contribution, k^^^nr^. As iV is 

' (JV-l) s™ 

increased, the scaling of this contribution with so 
becomes more and more rapid, aiding in the fit- 
ting of C2JV+2- The decrease in the coefficient fac- 
tor, 1/(A^ — 1), also means that the corresponding 
FESR is more strongly D — dominated, a de- 
sirable situation for the determination of ag . The 
latter effect is dominant for sufficiently large N. 
In addition, as N is increased, u'jv(y) ^ (1 ~ 
whose single zero at s = sq provides less strong 
suppression of contributions from the region of 
the timelike point on the OPE contour. We thus 
restrict our attention to the W2, • ■ ■ , FESRs. 

Since D = 2 contributions are negligible and 
D = A contributions are fixed by phenomenologi- 
cal input, the only OPE parameters to be fit us- 
ing the wn FESR are as(m^) and C2N+2- The 
C2N+2 will of course depend on the channel (V, A 
or V-l-A) being considered. The values for asim^) 
obtained using the different wn and/or different 
channels should, however, be consistent, and this 
consistency represents an important cross-check 
on the reliability of the analysis framework. 

We have analyzed the W2, • • ■ , ^6 FESRs using 
the final 2005 ALEPH isovector data and covari- 
ances, in each of the V, A and V-l-A channels. A 
similar analysis has been performed for the V+A 
channel using the OPAL data and covariances. 
See Ref. [5] for further details, and a discussion 
of the reasons for the analysis choices. 

We report here only on the results for as(m^). 
A full discussion of the errors, and results for the 
C2N+2, may be found in Ref. [5]. Results ob- 
tained using the CIPT prescription are presented 
in Table [H For each entry, the first error is exper- 
imental (computed using the experimental covari- 
ance matrix, and including the 0.32% normaliza- 
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Figure 1. Comparison of the Fy{so) corresponding to (i) our fits and (ii) the 2005 ALEPH fit, and 
various weights having degree < 3. The hght (heavy) dotted line corresponds to the ALEPH fit (our fit) 
for the hght (heavy) dashed hne to the ALEPH fit (our fit) for W2, the hght (heavy) dot-dashed 

line to the ALEPH fit (our fit) for ws, and the light (heavy) double-dot-dashed line to the ALEPH fit 
(our fit) for 'w{y) = y{l — y)^. 
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Table 1 

Results of the wn CIPT-based FESR fits for 
Q!s("t-t) obtained using either the ALEPH or 
OPAL data and covariances. The first error is 
experimental and the second theoretical. 



Data set 


Channel 


Weight 


as {ml) 


ALEPH 


V 


W2 


0.32f(7) 


(12) 






W3 


0.32f(7) 


(12) 






W4 


0.321(7) 


(12) 






W5 


0.321(7) 


(12) 






Wq 


0.321(7) 


(12) 




A 


W2 


0.319(6) 


(12) 






W3 


0.319(6) 


(12) 






Wi 


0.319(6) 


(12) 






W5 


0.319(6) 


(12) 






Wq 


0.319(6) 


(12) 




V+A 


W2 


0.320(5) 


(12) 






W3 


0.320(5) 


(12) 






Wi 


0.320(5) 


(12) 






W5 


0.320(5) 


(12) 






We 


0.320(5) 


(12) 


OPAL 


V+A 


W2 


0.322(7) 


(12) 






W3 


0.322(7) 


(12) 








0.322(7) 


(12) 






W5 


0.322(7) 


(12) 






We 


0.322(8) 


(12) 



tion uncertainty) , while the second is theoretical. 
The corresponding FOPT-based results may be 
found in Ref. [5]. 

From the table we see excellent consistency be- 
tween the ALEPH-based V, A and V+A results, 
as well as between the ALEPH and OPAL V+A 
results. There is also extremely good consistency 
within each channel between results obtained us- 
ing the different wn- This consistency is realized 
only after fitting the small, but non-negligible, 
D > 4 contributions on the OPE sides of the var- 
ious FESRs [5J. The FOPT results correspond- 
ing to different but the same channel display 
significantly less good consistency [5]. From the 
table, one sees that theoretical errors are a fac- 
tor of ~ 2 larger than the corresponding experi- 
mental errors. The dominant contribution to the 
theoretical error is that associated with the trun- 



cation of the series for the dominant D — OPE 
contribution. Our truncation error is the sum in 
quadrature of the shift in ag produced by includ- 
ing the last incorporated term in the truncated 
Adler function series and the difference between 
the results produced by the FOPT and CIPT 
evaluations. The latter contribution is the larger 
of the two and, as noted previously in the liter- 
ature, shows no signs of decreasing with increas- 
ing truncation order. The FOPT-CIPT difference 
thus dominates the current uncertainty. 

Our central determination for is based on 
the V+A analyses, which has the smallest ex- 
perimental errors. Averaging the ALEPH- and 
OPAL-based results using the non-normalization 
component of the experimental errors yields 

as{ml) = 0.3209(46)(118) (6) 

where the errors are experimental (including nor- 
malization) and theoretical, respectively. 

The Uf ^ 5 result, as{M^), is obtained from 
Eq. ((6|) using the standard self-consistent com- 
bination of 4-loop running with 3-loop matching 
at the flavor thresholds [15]. With mc{rnc) = 
1.286(13) GeV, rnb{mb) = 4.164(25) GeV [16], 
matching thresholds rmcbimcb) with r varying 
between 0.7 and 3, and standard estimates for 
the effect of the truncated running and matching, 
the evolution-induced uncertainty on as(M^) is 
0.0003. Our final result is then 

a,(Af|) = G.1187(3)e„oi(6)exp(15)th. (7) 

The difference between this value and that ob- 
tained in Ref. [H], 0.1212(11), serves to quantify 
the impact of the D > 8 contributions neglected 
in the earlier spectral weight analyses. 

The result, Eq. ([7]), is in excellent agreement 
with recent independent determinations, includ- 
ing that of the updated global EW analysis [H 
[TT and two recent updates [17|18j of the older 
HPQCD/UKQCD lattice determination [19]. 

With theory errors dominant, and the D = 
truncation uncertainty dominating the theory er- 
ror, further improvement will be possible only 
if a better understanding of the truncation un- 
certainty can be obtained. In a recent paper, 
Beneke and Jamin ^20J investigated this issue us- 
ing a model which incorporates the general struc- 
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ture associated with the first two IR renormalon 
and leading UV renormalon singularities of the 
resummed D = series. This particular version 
of what could be a more general model was used 
to argue in favor of the FOPT over the CIPT pre- 
scription. The FOPT evaluation, as well as the 
resummed model, together with assumed values 
for Cq and Cs, were also used to determine as- 
While it was shown in Ref. [5] that the assumed 
Cg and Cg values are not consistent with the ag 
obtained using FOPT, the underlying approach 
remains extremely interesting. In fact, one can 
see that the minimal version of the model em- 
ployed in Ref. [20] predicts CIPT approximations 
which deviate from the model version of the re- 
lated resummed series by amounts that are size- 
able, and a factor of '~ 2 larger for u''*'°-' than 
for W2 and W3 21] . The deviations are signifi- 
cantly smaller for the FOPT versions. Were the 
model to provide a good representation of the true 
resummed series, one would thus expect no set 
of as, Cg and Cs to provide a good quality si- 
multaneous CIPT fit for the W2, W3 and 
spectral integrals, while a reasonable quality si- 
multaneous fit would be expected to exist using 
the FOPT prescription. In fact, just the oppo- 
site occurs: the values of as, Cg and Cs ob- 
tained from a combined CIPT fit to the W2, W3 
FESRs provide also an excellent representation of 
the w'--^^^ spectral integrals, while those obtained 
from the corresponding FOPT fit provide a very 
poor one. This suggests an extended version of 
the Beneke-Jamin model, perhaps taking into ac- 
count more IR and/or UV renormalon singular- 
ities, is likely to be needed. The fact that it is 
possible to reach such a conclusion, however, im- 
mediately makes clear that the sets of variously 
weighted so-dependent spectral integrals provide 
significant constraints for use in constructing such 
generalizations. It is thus likely that such mod- 
elling can be improved in future, and used to re- 
duce the truncation uncertainty component in the 
determination of as ■ 
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